Electron-boson glue function derived from electronic Raman scattering 
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Raman scattering cross sections depend on photon polarization. In the cuprates nodal and antin- 
odal directions are weighted more strongly in Big and Bi g symmetry, respectively. On the other 
hand in angle-resolved photoemission spectroscopy (ARPES), electronic properties are measured 
along well-defined directions in momentum space rather than their weighted averages. In contrast, 
the optical conductivity involves a momentum average over the entire Brillouin zone. Newly mea- 
sured Raman response data on high-quality Bi2Sr2CaCu2 0s+i single crystals up to high energies 
have been inverted using a modified maximum entropy inversion technique to extract from Bi g 
and B2 g Raman data corresponding electron-boson spectral densities (glue) are compared to the 
results obtained with known ARPES and optical inversions. We find that the Big spectrum agrees 
qualitatively with nodal direction ARPES while the B\ g looks more like the optical spectrum. A 
large peak around 30 — 40meV in Bi 9 , much less prominent in -829, is taken as support for the 
importance of (tt, it) scattering at this frequency. 

PACS numbers: 74.20. Mn, 74.72.Gh, 78.30.-j 



Boson structures seen in tunneling and to a lesser ex- 
tent in optics in conventional superconductors have given 
us detailed information about the electron-phonon inter- 
action in these materials. While similar structures have 
been identified in the cuprates in tunneling, point contact 
junctions 1 as well as scanning tunneling spectroscopy 
(STS) 2 ' 3 , in angle-resolved photo emission (ARPES) 4 " 10 , 
and, particularly, in optics 11-14 , some of the details as- 
sociated with the recovered electron-boson spectral den- 
sity remain controversial, particularly the nature of the 
bosons involved. Some investigators stress the role of 
phonons 1-7 while others favour spin fluctuations 8-14 . Of 
course, both mechanisms are expected to contribute to 
the effective electron-boson interaction so that the de- 
bate really centers on which might be dominant 14-16 . 
Recent ARPES data have found a 6% softening of the 
boson 'kink' in the nodal direction renormalized elec- 
tronic dispersion curves for Bi 2 Sr2CaCu208+5 (Bi2212) 
compounds on substitution of 16 by 18 15 ' 16 . Based 
on an examination of the detailed shape of the electron- 
boson spectral function, 7 2 x(u;), (glue function) obtained 
from a maximum entropy inversion 17 of data on a simi- 
lar sample 10 of Bi2212, Schachinger et al. 15 have argued 
that the isotope substitution data can be understood if a 
peak seen in the electron-boson spectral density I 2 \{uj) 
around 65meV which contains about 10% of the total 
area is assigned to an oxygen phonon mode. This leaves 
90% of the effective spectrum which extends to 400 meV, 
much larger than any phonon energy, to possibly come 
from spin fluctuations. 

So far / 2 x(w) has been recovered from tunneling 1,18 , 
nodal direction ARPES 6-8,10 , and optical data in the 
cuprates 11-14 ' 19,20 . In principle, I 2 x(w) is anisotropic 



and it will be different for each momentum direction. 
Also quasiparticle and optical spectral densities will not 
be the same. While both involve the same bosons differ- 
ent weighting electronic factors apply. This is also true 
for Raman for which different vertices apply for different 
photon polarization 21-26 . Nevertheless, it is important 
to obtain the corresponding electron-boson spectral den- 
sity associated with B lg and B 2g Raman data and to 
understand how these might differ from those obtained 
from ARPES and from optics and to establish points of 
consistency between these various spectra. Based on a 
memory function approach to the Raman cross section 25 
one can extract a corresponding Raman scattering rate 
which depends on polarization. In this work we apply 
a maximum entropy inversion technique to extract from 
such data a spectral density in analogy to what is done 
for the infrared conductivity. 

The analysis starts from new experimental results 
from freshly prepared high-quality single crystals of 
Bi2Sr2CaCu208+,5. The crystals were grown in a mirror 
furnace using the traveling-solvent floating-zone (TSFZ) 
technique and post-annealed at 870 °C to arrive at opti- 
mal doping (p = 0.16 holes/Cu0 2 ) with a T c of 94.5 K 
and a transition width below 1 K. The quality of the sam- 
ples is crucial since defects lead to a strongly enhanced 
cross section at high energies 27 . 

The basis of the analysis is the Raman response func- 
tion R^xZiTfCj), where \x represents the scattering sym- 
metry. Here we focus on /1 = B\ g and £>2g- The re- 
sponse is derived from the measured cross section as 
R^xl = {io I /ios)[l + n^,T)]- 1 d 2 aJ{d^dLUs) 26 . Here, 

(S) is t ne frequency of the incident (scattered) light. 
f2 is the solid angle of acceptance of the collection optics 
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and n(ui,T) — {exp[hu>/ (fcgT) — l]} -1 is the Bose factor. 
<7^ is corrected for the sensitivity of the instrument. Since 
the optical constants vary only little for visible light, in- 
terface effects can be absorbed in the constant the 
magnitude of which is irrelevant here. For revealing the 
pure symmetry components of the response in a reliable 
fashion we measured spectra at all six main polarizations 
of the afe-plane. This allows us to check the internal con- 
sistency of the procedure. Thus, the B\ g response, for 
example, reads in a short hand notation 



B 



is 



xx + xy + RL-^ (RR + xy + xx 



(1) 



where x = [100], x' = [110], R = [100] + i[010], etc. The 
typical number of photon counts per point is between 
1000 and 10000 which results in a relative statistical error 
in the range of one to three percent. For the experiments 
shown here the read-out noise of the CCD detector is 
negligable. 

Opel et al. 25 show how a Raman scattering 
rate, r M (T, cj), can be extracted quantitatively from 
i?^X^(T, (jj). At a given symmetry the sensitivity in mo- 
mentum space is uniquely defined, and the main con- 
tributions in B\g and B 2g symmetry come from the re- 
gion close to (it, 0) (and equivalent points of the Bril- 
louin zone) and from the center of the quadrant, re- 
spectively. The derived scattering rates r M (T, w) are 
analogous to the optical scattering rate, r~ pt (T,w), ob- 
tained from infrared measurements 11 but correspond to 
different parts of the BZ dictated by symmetry /i 26 . 
In the normal state and within a Kubo formalism we 
can show that the Raman spectrum 24 and the opti- 
cal response 28-32 is related to the appropriate electron- 
boson spectral density, 7 2 x(^), to a good approximation 
through the equation 33-35 

oo 

IV(:Z» - r~} mv = J dvK{uj,v;T)ll X {v), (2) 
o 

where r~l mp is a Raman impurity scattering rate and 
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where T denotes the temperature. For the conductiv- 
ity which involves an average over all momentum direc- 
tions 9 the appropriate electron-boson spectral density 

is l lx{ v ) = ktfx{ v i$))o with (•••>0 thc average over 
the directions 9. For the Raman case with symmetry 
H there is an additional weighting of cos 2 (26*) for Bi g 
and sin 2 (20) for B 2g - Here, I^xi 11 , 9) is the electron-boson 
spectral density associated with momentum direction 9. 
There is another difference between quasiparticle and 
transport quantities associated with vertex corrections. 



These are expected to mainly change the magnitude of 
the distribution functions with shape changes secondary. 

Jiang and Carbotte 24 give the formula for the lowest 
order Raman susceptibility for an interacting electron 
system in the form 

- T ^2^2 Tr [7£(q)T30(q,M*>m)T3<?(q,iu> m - iu n )] , 



m q 



(4) 



with uj m and v n the fcrmionic and bosonic Matsub- 
ara frequencies, respectively, q the momentum vector, 
T3 the third Pauli matrix and Q(q,iuj m ) the electronic 
matrix Green's function in Nambu notation. Equa- 
tion (4) is valid in the superconducting as well as nor- 
mal state. Here we consider only the latter. Equa- 
tion (4) differs from the well known formula for the op- 
tical conductivity 24 ' 28-31 only through the factor 7 2 (q), 
Ref. 24 , which is to be replaced by e 2 Vp x (q) where VF,x{q) 
is the x-component of the Fermi velocity at momentum q 
which in the free electron model is assumed constant and 
Vp x (q) = Vp/2 in two dimensions. For infinite bands 
with constant electronic density of states N(0), Eq. (4) 
can be written in a more convenient form as 
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This equation has the same form as Eq. (22) of Sharapov 
and Carbotte 35 reported in their study of the effect of the 
energy dependence of the quasiparticle density of states 
on the far-infrared absorption in underdoped cuprates. 
Equation (5) differs as there is now an integration over 
angles 9, the additional Raman vertex 7 2 (9), and a differ- 
ent numerical factor. The same algebraic manipulation 
as is applied in Ref. 33 gives: 



<w_ 

2tt 



de[f(e)-f(e + co)] 



7 *(0)Im [£„(e + w,0)-E*(e,0)] ,(6) 



where (e, 9) is the quasiparticle self energy due to the 
directional electron-boson spectral density J 2 x(^, 0) and 
the * indicates the complex conjugate. This leads directly 
to our fundamental Eq. (2) when we use the relationship 
for the imaginary part of the quasiparticle self energy in 
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terms of I 2 x(Q, 6), namely 



ImE M (w,0) 



<milx{si,e) 



2coth 



. uj + . 
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In the approximation of Eq. (5) for x M (w) we can con- 
struct the corresponding Raman cross section from 



(8) 



where ujX^uj) is the Kramers-Kronig transform (KK) of 
r AJ (w) as described by Opel et al. 25 [Eq. (A8)] to whom 
we refer for a detailed discussion of how the scattering 
rate r Al (w) is extracted from the data on the Raman 
cross section corresponding to B\ g and B 2g symmetries. 
In Fig. 1 we present our results for the electron-boson 
spectral density / 2 x(w), Fig- l( a ) f° r Big (antinodal) 
and Fig. 1(b) B 2g (nodal) polarizations. We show for 
comparison equivalent results obtained previously from 
optics 17 ' 36 , Fig. 1(c), and from nodal direction ARPES, 
Ref. 10 , Fig. 1(d). In Fig. 2(a) we show our maximum 
entropy fits to the B 2g Raman scattering rates for two 
temperatures. The light solid (blue) dots are experiment 
at T = 98 if and the heavy solid (black) curve our fit. 
The light open (black) dots give the data for T = 298 K 
with the heavy dashed (red) curve our fit. Note that the 
low frequency part of the data at 98 K is almost linear, 
a feature which is well captured by the theoretical fit 
except for the very lowest u> region. On the whole the 
general trend is well described by our data reconstruc- 
tion. To compensate for the rather big zero frequency 
offset of TB 2g ( w - ► 0, T) an impurity scattering rate of 
Tn 1 . „ = 44meV had to be introduced in Eq. (2). In 
Fig. 2(b) we show a comparison between data and theo- 
retical fit for the B 2g Raman cross section (see Ref. 25 for 
notation) which is the measured quantity from which the 
Raman scattering rate Tb 2!I (T,uj) of Fig. 2(a) has been 
extracted. To get this quantity we applied Eq. (8), with 
As 2g (T, u>) obtained from a Kramers-Kronig transform 
of our fitted Raman scattering rate Tb 2!I {T,uj). We see 
good agreement except for a low frequency peak in the 
theoretical curve at T — 98 K [heavy solid (black) curve] 
not present in the data. This peak is traced to our use of 
an infinite band approximation. The experimental data 
on r M (T, uj) develop peaks around 700 me V (not shown 
here) for both B lg and B 2g symmetry which is under- 
stood to be the result of a reduction in the electronic 
density of states as a band edge 10 is approached. This 
effect is not captured in our calculations and translates 
into important differences between calculated and exper- 
imental A M (T, uj). Furthermore, we observe that theory 
deviates consistently from experiment to lower values for 
energies greater that ~ 350 meV. This can also be traced 
to our use of the infinite band width approximation. A 



more complete discussion of the effect of finite bands in 
the inverted electron-boson spectrum can be found in 
Ref. 10 . 

Similar results are found for the B\ g data reconstruc- 
tion. Vb 1 (uj,T) is reproduced equally well except for 
some deviations at low energies seen in the T = 302 K 
data because of a linear frequency dependence in the ex- 
perimental TB lg (co,T) data extending from <~ 100 meV 
down to uj = which cannot be reproduced by theory. 
Such a low energy linear dependence at room tempera- 
ture has not been observed in the B 2g data [light open 
(black) dots in Fig. 2(a)]. The zero frequency offset of 
Tb 1s {oj — > 0, T) was compensated by an impurity scat- 
tering rate r^ 1 imp = 81.6 meV about twice as much as 
was necessary for B 2g . The data reconstruction of the 
Bi g Raman cross section also reveals deviations from ex- 
periment at low and high frequencies which, again, can 
be understood to be the result of our infinite band ap- 
proximation on which Eq. (2) is based. 

We turn now to a comparison of antinodal with nodal 
Raman results of Figs. 1(a) and 1(b). The shapes of the 
distributions obtained are quite distinct. For B\ g there is 
a large peak at ~ 29 meV in the T = 98 K spectrum fol- 
lowed by a dip and then a second peak around 300 meV. 
As the temperature is increased there is a clear evolution 
of the spectrum with the low energy peak decaying in 
amplitude, broadening and moving towards higher ener- 
gies. The same trend, although less pronounced, is ob- 
served for the second peak. The valley between the peaks 
becomes progressively filled in but it still very much re- 
mains, even for T = 302 K although the effects are much 
less pronounced. This is also seen when the B 2g (nodal) 
Raman spectrum is considered instead of B\ g (antin- 
odal). The shape of I 2 x(w) for B\ g as well as its change 
with increasing temperature agrees well with previous 
trends for the optical case. In Fig. 1(c) we show results at 
T = 100 K, 200 K and 295 K in a sample of Bi2212 based 
on data by Tu et al. 36 . The prominent peak is at a posi- 
tion (~ 44 meV) slightly different from our Raman spec- 
trum, but the overall shape at 100 K is in good qualitative 
agreement with the Raman result [Fig. 1(a)] including the 
second peak, the valley between the two peaks, and the 
temperature evolution. A more extensive set of data on 
a similar sample but doped with some yttrium is found 
in Ref. 12 . The data presented in their Fig. 2, top frame, 
is also in good qualitative agreement with our Raman 
results. The important observation is that both antin- 
odal Raman and optical data show a prominent peak in 
the electron-boson spectral density around 30 to 40 meV 
which does not appear in ARPES and appears much less 
prominently in nodal Raman data. More specifically, the 
fractional area under this peak is only ~ 3% in B 2g as 
compared with ~ 23% in B\ g . As noted above, these 
two structures evolve with temperature in the same way 
which indicates their common origin. Part of this tem- 
perature evolution could be due to the reduced resolution 
intrinsic to our unbiased inversion method as the temper- 
ature is increased. (See Ref. 17 .) This is consistent with 
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FIG. 1: (Colour on-line) The electron-boson spectral density (dimensionless) as a function of energy ui in meV from Bi g (a) 
and B2 g (b) Raman data. We show in (c) results 17 obtained from optical data 36 and in (d) a result 10 obtained from nodal 
direction ARPES 6 . 
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FIG. 2: (Colour on-line) (a) Fit to the E>2 g Raman data for the scattering rate F(cj,T) from our maximum entropy reconstruc- 
tion, (b) The corresponding fits to the measured £?2 9 cross section. 



the results obtained with the biased inversion method 
discussed later in Fig. 3 which shows less temperature 
dependence. 

For nodal ARPES there is no peak at ~ 30meV but 
one appears instead at much higher energy ~ 65 meV. 
This higher energy has often been identified with cou- 
pling to oxygen phonons 4-6 ' 15,16 and we will return to 
this issue later. Such a peak is not seen in our Raman 
spectra which indicates that such effects are mainly con- 
fined to the nodal direction and do not appear in averaged 
quantities such as B\ g and even B^ g Raman even though 
for this latter polarization the Raman vertex peaks in the 
nodal direction. The fact that the peak at ~ 30meV is 
stronger in B\ g (antinodal) than in B2 g (nodal) and is not 



seen in nodal ARPES is consistent with a boson which 
is associated with scattering through momentum trans- 
fer q of (ir,ir). Such a vector corresponds to transitions 
between those parts on the Fermi surface which lie also 
on the antiferromagnctic Brillouin zone and is, therefore, 
closer to the antinodal direction. Consequently, quanti- 
ties that emphasize the area around (7r, 0) such as B\ g 
Raman and, to some extent, optics should show a strong 
peak. In contrast, the peak at ~ 30meV is expected to 
be weaker in B2 g symmetry weighing out the nodal part. 
This is what we observe and have shown in the data of 
Figs. 1(a) and (b). 

Instead of using a maximum entropy technique to 
invert Eq. (2) van Heumen et a/. 14 used a histogram 
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to characterize the electron-boson spectral density, 
I^ pt x(uj), derived from optics. The histogram is then used 
to directly reconstruct the experimental data by inversion 
of Eq. (5) modified for the optical conductivity. They 
find less temperature dependence than we have here, for 
the position and width of the peak around 30meV in 
Bi g . So far we used an unbiased maximum entropy in- 
version of Eq. (2) in which the default model 17 is set to 
a constant at all temperatures. Another method is the 
so-called biased maximum entropy inversion in which the 
default model is set to the previous next lower tempera- 
ture solution. We can expect the solutions of these two 
methods to be different because the inversion of Eq. (2) 
is an ill-posed problem which usually has more than one 
solution. It was pointed out by Yang et al. 19 that in the 
case of HgBa 2 CuC>4 +( 5 the biased inversion of optical data 
resulted in electron-boson spectra very similar to those 
reported by van Heumen et al. 14 for the same material. 

Results of such a biased inversion of the Raman B\ g 
and Big data are shown in Figs. 3(a) and 3(b), respec- 
tively. We see that the peak at ~ 30meV for the B\ g 
spectra shifts less with temperature than is shown in 
Fig. 1(a) although it still loses amplitude with increas- 
ing temperature T. This points to the possibility that at 
least part of this peak comes from a phonon, and its con- 
tribution could be determined from the amplitude of this 
peak at room temperature. The situation is quite differ- 
ent for the B 2g polarization. Comparison of the spectra 
shown in Fig. 3(b) with the ones presented in Fig. 1(b) 
reveals that there is very little difference between the two 
sets of spectra. In particular, the position of the low en- 
ergy peak (~ 40meV at 98 K) shows almost the same 
temperature dependence and is almost smeared out at 
room temperature. Furthermore, the inversion method 
has little or no influence on the shape and size of the 
high energy part (lj > 100 meV) of the B lg and B 2g spec- 
tra. Figures 3(c) and 3(d) show results for the often used 
mass enhancement parameter A M = 2 du> I^x{uj) / 'ui as 
a function of temperature for B lg and B 2g polarization, 
respectively. There is very little difference between the 
two methods of inversion. 

Returning to Fig. 1(b), the nodal direction Raman 
spectrum (B 2g ) has a very different characteristic shape 
as compared to B\ gi Fig. 1(a). While it can be character- 
ized as also having two peaks, the valley between them is 
not pronounced and its spectral weight is much more uni- 
formly distributed below 300 meV with a relatively sharp 
drop off beyond this energy. This shape is much closer to 
what has been found in inversions of the nodal direction 
ARPES data reproduced in Fig. 1(d). This is expected 
since the B 2g Raman vertex peaks at the nodal direction. 
Note that ARPES is strictly directional and samples only 
the nodal direction while B 2g Raman probes an extended 
part of the BZ weighted by sin 2 (26>). Nevertheless, the 
agreement as to shape between B 2g Raman and nodal di- 
rection ARPES gives one confidence that both methods 
are measuring the same boson spectrum. Note that there 
is nothing which limits the application of this method to 



the high-T c cuprates and our inversion technique has a 
more general applicability to other metals. 

While our inversions provide us with a good handle 
on the size and qualitative shape of the spectral density, 
the question as to the origin of the boson involved is 
more difficult to answer in a definite way. Certainly if 
phonons are involved we would expect 7^x( a; ) to mirror 
the phonon frequency distribution while for spin fluctua- 
tions we should see an image of the local spin susceptibil- 
ity. The fact that J 2 x( w ) shows very significant spectral 
weight up to 400 meV means that excitations other than 
phonons having energies below 100 meV are involved. For 
spin fluctuations the energy scale is set by the exchange 
coupling, J, which enters, e.g., the t — J model as a pa- 
rameter, and this is consistent with the large energy scale 
seen here. It is also consistent with recent numerical 
studies of the t — J model by Maier et al. 37 in which an 
effective electron-boson spectral density associated with 
short range spin fluctuations is extracted and identified 
as the pairing glue. It displays many of the features seen 
in our empirical spectra. The cellular dynamical mean- 
field calculations of Kyung et al. 38 based on the Hubbard 
model also give qualitatively similar results for the spec- 
tral density and offer further microscopic support for an 
interpretation of our derived spectra as due largely to 
spin fluctuations. 

On the other hand one does expect, and experimen- 
tal data provide support for some contribution from the 
electron-phonon interaction 3 ' 5-7 . While the change in 
critical temperature T c on substitution of 16 — > 18 O 
is small for optimally doped samples, it is nonzero and 
it can be large for the underdoped case. However, the 
latter fact can also be understood as due to an energy 
dependence in the electronic density of states 39 or to a 
pseudogap formation 40 while at the same time the under- 
lying contribution of phonons to the pairing interaction 
remains small. Recent ARPES experiments along the 
nodal direction 16 found a shift in the 'boson kink' in the 
renormalized dispersion curves of Bi2212 upon oxygen 
isotope substitution. This was assigned by Schachinger et 
al. 15 as a 10% phonon contribution to the electron-boson 
spectral density found in inversions 10 of ARPES data as 
reproduced in Fig. 1(d). They assign the peak around 
65meV to electron-phonon effects and the rest due to 
spin fluctuations. Accounting for finite band effects this 
translates into a contribution to the quasiparticle mass 
enhancement due to phonons of A ~ 0.2 at a temperature 
of T = 17 K. A similar estimate was found by Devereaux 
et al. 41 due to the buckling and breathing phonon modes. 
Recent calculations based on local-density approximation 
also predict mass enhancements due to phonons 42-45 but 
are up to an order of magnitude smaller than our result. 

The aim of this work was to extract from the electronic 
Raman cross section in B\ g and B 2g (/i) symmetry, the 
corresponding electron-boson spectral density I 2 x(u;). 
Comparison with equivalent, previously extracted forms 
from optical and nodal ARPES data shows remarkable 
consistency between the results obtained with such dif- 



6 




FIG. 3: (Colour on-line) (a) The electron-boson spectral density I 2 x(u) (dimensionless) as a function of energy u in meV from 
Big Raman scattering rates using a biased maximum entropy inversion as described in the text, (b) The same as (a) but now 
for B29 symmetry. Frames (c) and (d) compare the temperature dependences of the mass enhancement parameters A obtained 
from biased and unbiased inversion of the B\ g and i?2g Raman scattering rates, respectively. 



ferent probes. Comparison of these results provides in- 
formation on the angular variation of the electron-boson 
self energy and corresponding spectral density around the 
Fermi surface. This arises because optics involves an av- 
erage over the entire Fermi surface while Raman weights 
predominantly the antinodal and the nodal directions for 
Big and B2 g symmetry, respectively. Only ARPES is 
perfectly directional. While we now have a good handle 
on the size and qualitative shape of the spectral den- 
sity, the question of what bosons might be involved is 
more difficult. Certainly, the general shape of the 7 2 x(w) 
should reflect the phonon frequency distribution if the 
electron-phonon interaction is dominant while it should 
mirror the shape of the local spin susceptibility if it is, in- 
stead, the spin fluctuations that are dominant. The fact 
that all spectral functions obtained have very significant 
spectral weight beyond 100 meV makes it clear that ex- 



citations other than phonons are primarily involved. An 
energy scale set by J is more consistent with our results. 
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